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Abstract. The discrete Fourier analysis on the 30°-60°-90° triangle is deduced from the 
corresponding results on the regular hexagon by considering functions invariant under the 
group G2, which leads to the definition of four families generalized Chebyshev polynomials. 
The study of these polynomials leads to a Sturm-Liouville eigenvalue problem that contains 
two parameters, whose solutions are analogues of the Jacobi polynomials. Under a concept 
of m-degree and by introducing a new ordering among monomials, these polynomials are 
shown to share properties of the ordinary orthogonal polynomials. In particular, their 
common zeros generate cubature rules of Gauss type. 
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1 Introduction 

In our recent works [9, 10, 11] we studied discrete Fourier analysis associated with translation 
lattices. In the case of two dimension, our results include discrete Fourier analysis of expo- 
nential functions on the regular hexagon and, by restricting to symmetric and antisymmetric 
exponentials on the hexagon under the reflection group A2 (the group of symmetry of the regu- 
lar hexagon), the generalized cosine and sine functions on the equilateral triangle, which can 
also be transformed into the generalized Chebyshev polynomials on a domain bounded by the 
hypocycloid. These polynomials possess maximal number of common zeros, which implies the 
existence of Gaussian cubature rules, a rarity that is only the second example ever found. The 
first example of Gaussian cubature rules is connected with the trigonometric functions on the 
45°-45°-90° triangle. The richness of these results prompts us to look into similar results on 
the 30° -60° -90° triangle in the present work. This case is also considered recently in [13] as 
an example under a general framework of cubature rules and orthogonal polynomials for the 
compact simple Lie groups, for which the group is Gi- 

It turns out that much of the discrete Fourier analysis on the 30°-60°-90° triangle can be 
obtained, perhaps not surprisingly, though symmetry from our results on the hexagonal domain. 
The most direct way of deduction, however, is not through our results on the equilateral triangle. 
The reason lies in the underline group G2, which is a composition of A2 and its dual A%, the 
symmetric group of the regular hexagon and its rotation. Our framework of discrete Fourier 
analysis incorporates two lattices, one determines the domain and the other determines the space 
of exponentials. Our results on the equilateral triangle are obtained from the situation when 
both lattices are taken to be the same hexagonal lattices [9] . Another choice is to take one lattice 
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as the hexagonal lattice and the other as the rotation of the same lattice by 90° degree [10], with 
the symmetric groups A 2 and A%, respectively. As we shall see, it is from this set up that our 
results on the 30°-60°-90° triangle can be deduced directly via symmetry. The results include 
cubature rules and orthogonal trigonometric functions that are analogues of cosine and sine 
functions. There are four families of such functions and they have also been studied recently 
in [13, 18]. While the results in these two papers concern mainly with orthogonal polynomials, 
our emphasis is on the discrete Fourier analysis and cubature rules, and on the connection to 
the results in the hexagonal domain. 

The generalized cosine and sine functions on the 30°-60°-90° triangle are also eigenfunctions 
of the Laplace operator with suitable boundary conditions. There are four families of such 
functions. Under proper change of variables, they become orthogonal polynomials on a domain 
bounded by two curves. However, unlike the equilateral triangle, these polynomials do not form 
a complete orthogonal basis in the usual sense of total order of monomials. To understand 
the structure of these polynomials, we consider the Sturm-Liouville problem for a general pair 
of parameters a, f3, with the four families that correspond to the generalized cosine and sine 
functions as a = +i, (5 = ±|. The differential operator of this eigenvalue problem has the form 

£ a ,/9 := -An(x,y)dl - 2A 12 (x,y)d x d y - A 22 {x,y)d 2 y + B 1 (x,y)d x + B 2 {x,y)d y . 

Such operators have long been studied in association with orthogonal polynomials in two va- 
riables; see for example [6, 7, 8, 16], as well as [1] and the references therein. Our operator £ a> p, 
however, is different in the sense that the coefficient functions Aij are usually assumed to be 
of quadratic polynomials to ensure that the operator has n + 1 polynomials of degree n as 
eigenfunctions, whereas A 2y2 hi our L a ^ is a polynomial of degree 3 for which it is no longer 
obvious that a full set of eigenfunctions exists. Nevertheless, we shall prove that the eigenvalue 
problem C a pu = Xu has a complete set of polynomial solutions, which are also orthogonal 
polynomials, analogue of the Jacobi polynomials. Upon introducing a new ordering among 
monomials, these polynomials can be shown to be uniquely determined by their highest term 
in the new ordering. As a matter of fact, this ordering defines the region of influence and 
dependence in the polynomial space for each solution. Furthermore, it preserves the m-degree 
of polynomials, a concept introduced in [13], rather than the total degree. In the case of a = +^ 
and /3 = +2, the common zeros of these polynomials determine the Gauss, Gauss-Lobatto and 
Gauss-Radau cubature rules, respectively, all in the sense of m-degree. It is known that the 
cubature rule of degree 2n — 1 exists if and only if its nodes form a variety of an ideal generated 
by certain orthogonal polynomials. It is somewhat surprising that this relation is preserved 
when the m-degree is used in place of the ordinary degree. 

The paper is organized as follows. The following section contains what we need from the 
discrete Fourier analysis on the hexagonal domain. The results on the 30°-60°-90° triangle is 
developed in Section 3, which are translated into generalized Chebyshev polynomials in Section 4. 
The Sturm-Liouville problem is defined and studied in Section 5 and the cubature rules are 
presented in Section 6. 

2 Discrete Fourier analysis on hexagonal domain 

Before stating the results on the hexagonal domain, we give a short narrative of the neces- 
sary background on the discrete Fourial analysis with lattice as developed in [9, 11]. We refer 
to [2, 3, 12, 14] for some applications of discrete Fourier analysis in several variables. 

A lattice L in R d is a discrete subgroup L = La '■= A7j d , where A, called a generator matrix, 
is nonsingular. A bounded set Q of called the fundamental domain of L, is said to tile R 
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with the lattice L if Q + L = ~R d , that is, 

^ Xn(x + a) = 1, for almost all x e M. d , 

aeL 

where xn denotes the characteristic function of Q. For a given lattice La, the dual lattice is 
given by L\ = A~ tr 1, d . A result of Fuglede [5] states that a bounded open set Q, tiles M. d with 
the lattice L if, and only if, {e 2ma ' x : a e L^} is an orthonormal basis with respect to the inner 
product 

(f,g)n = ^ T ^f(x)g{x)dx. (2.1) 

Since L\ = A- tr Z d , we can write a = A~ tr k for a e L\ and k e Z d , so that e 2 ™' x = e ^ ikt ' A ' lx . 

For our discrete Fourier analysis, the boundary of VL matters. We shall fix an VL such that 
e f2 and ft + AL d = ~R d holds pointwisely and without overlapping. 

Definition 2.1. Let FIa and be the fundamental domains of A7L d and BIj d , respectively. 
Assume all entries of the matrix N := B tr A are integers. Define 

A N := {keZ d : B- tr keQ A } and A^ := {k e Z d : A~ tr k e fl B }. 

Furthermore, define the finite-dimensional subspace of exponential functions 

V N := span{e 2 ^ fctrA ~ lx , k e A^}. 

A function / defined on R d is called a periodic function with respect to the lattice A7L d if 

f(x + Ak) = f(x) for all k e Z d . 

The function x 1— > e 2mkt ' A lx is periodic with respect to the lattice AX d and Vat is a space 
of periodic exponential functions. We can now state the central result in the discrete Fourier 
analysis. 

Theorem 2.2. Let A, B and N be as in Definition 2.1. Define 

for f, g in C{VLa), the space of continuous functions on VLa- Then 

(f,g)n A =(f,g) N , f,geV N . (2.2) 

It follows readily that (2.2) gives a cubature formula exact for functions in Vn- Furthermore, 
it implies an explicit Lagrange interpolation by exponential functions, which we shall not state 
since it will not be needed in the present work. 

In the following, we shall call the lattice La as the lattice for the physical space, as it 
determines the domain on which our analysis lies, and the lattice Lb as the lattice for the 
frequency space, as it determines the points that defines the inner product. 

The classical discrete Fourier analysis of two variables is the tensor product of the results 
in one variable, which corresponds to A = B = J, the identity matrix. We are interested in 
choosing A as the generating matrix H of the hexagonal domain, 

H = (-1 2) With M H = {xeM 2 : -l^x 2 ,^f±±f <1 
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If we choose B = ^H, so that N = B tr A has all integer entries, we are back to the situation 
studied in [9], which is the one that leads to the discrete Fourier analysis on the equilateral 
triangle. The other choices are considered in [10]. 

For the case that we are interested in, we choose A = H, the matrix for the hexagonal lattice 
in the physical space, and B = nH~ tr with neZ, the matrix for the hexagonal lattice in the 
frequency space. Then N = B tr A = nl has all integer entries. This case was studied in [10], 
which will be used to deduce the case that we are interested in by an additional symmetry. As 
shown in [9, 17], it is more convenient to use homogeneous coordinates (ti,t 2 ,t 3 ) defined by 



/ V3 
2 



\ _vi 

\ 2 




Ex, 



(2.3) 



which satisfy £1 + ^2 + ^3 = 0. We adopt the convention of using bold letters, such as t to denote 
points in homogeneous coordinates. We define by 



l 3 H ■= {t = (ti,t 2 ,h) eM 3 :ti + t 2 + t 3 = 0} and H T := I? 



n R 3 H 



the spaces of points and integers in homogeneous coordinates, respectively. In such coordinates, 
the fundamental domains of the lattices La and L b are then given by 



n : = n A = {t e : 

Q B = {t e 
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-l<h,t 2 ,-t 3 ^ 1}, 

-n < h - t 2 ,h - t 3 ,t 2 - t 3 ^ 11} , 



where Qa can be viewed as the intersection of the plane t\ + 1 2 + t 3 = with the cube [— 1, l] 3 . 
Define the index sets in homogeneous coordinates 

M n := {j e rf : —n j u j 2 , j 3 < n, j = (mod 3)}, 
El := {ketf : —n ^ k 3 -k 2 ,h - k 3 ,k 2 -k t ^ n}, 

where t = (mod m) means, by definition, t\ = t 2 = t 3 (mod m). We note that H n and ll 
serve as the symmetric counterparts of A^r and JV N , respectively, so that M n determines the 
points in the discrete inner product and Mn determines the space of exponentials. Moreover, 



the index set 



can be obtained from a rotation of Mn, as shown in the following proposition. 
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(-i,o,i)< 




(1,0,-1) 



(0,-1,1) (1,-1,0) 

Figure 2.1. Q,a in Cartesian coordinates (left) and homogeneous coordinates (right). 



Proposition 2.3 ([10]). For t = (h,t 2 ,t 3 ) e R 3 H , define t := (t 3 - t 2 ,t x - t 3 ,t 2 - h). Then 



irt 



i/keH n and k e H n i/ke 
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Figure 2.2. fi^ in Cartesian coordinates (left) and homogeneous coordinates (right). 

Proposition 2.3 states that HI n = : = {k : keil}. Similarly, we can define EI = : = 
{k : k e H^} = {j e tf : j = (mod 3)}. The set is the index set for the space of 
exponentials. Define the finite-dimensional space Tin of exponential functions 

%l : = span j<£ k = e^* :kei], 

By induction, it is not difficult to verify that 

dirndl = 



|H n | 



n 2 + n + 1, if n ^ 1 (mod 3), 
n? + n — 1, if n = 1 (mod 3). 



Under the homogeneous coordinates (2.3), x = y (mod H) becomes t = s (mod 3). We call 



(-n,n,0) (0,n,-n) 



(-n,n,0) (0,n,-n) 



(-n,0,n)^ 



(0,-n,n) 




(n, 0,-n) (-n,0,n)< 



i,-n,0) 



(-n,re,0) (0,n,-n) 



)(n,0,-ra) (-ii,0,n 




>(n,0,-re) 



(0,-ii,ii) (n,-n,0) 



(0,-n, n) (n,-n,0) 



Figure 2.3. H„ for n = 9 (left), n = 10 (center) and n = 11 (right). 



(-o,2o,-o) 



(-2a, a, a 



(-a,-a,2a)l. 




-a, 2a, -a) 



(2a, -a, -a) (-a, -a, 2a) 




(-a,2a,-a) 



(a,a,-2a) (-2a,a,a)f-; >i(a,a,-2a) (-2a, a, a) 



(2a, -a, -a) (-a, -a, 2a) 



(a, -2a, a) 



(a. -2a, a) 




(a,a,-2a) 



(2a, -a, -a) 



a, -2a, a) 



Figure 2.4. Hj, for n = 9 (left), n = 10 (center) and n = 11 (right), where 



a function / H-periodic if /(t) = f(t + j) whenever j = 0(mod3). Since j,k e H implies that 
2 j • k = (ji — j2)(^i — ^2) + 3J3/C3, we see that (/>j is H-periodic. 
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Theorem 2.4 ([10]). The following cubature rule holds for any f e T-i\ n _i, 
i r i f 1 ' jeH - 

ml /(*)* = ^E4 B) /(i). 4 n) = i j^n, (2.4) 



J 6 



where H° , and denote the set of points in interior, set of vertices, and set of points 
on the edges but not on the vertices; more precisely, M° n = { j e H : — n < ji, j'2, J3 < n}, = 
{(n, 0, -n)a e M : a e A 2 } and W n = M n \(W n u B%) = {(j, n - j, -n)a e M : 1 < j < n - 1}. In 
particular, let Q n f denote the right hand side of (2.4); then for any k e rf, Q n 0k = 1 i/k = 
(mod 3n) and Q n 4>k = otherwise. 

Here we state the main result in terms of the cubature rule (2.4), from which the discrete 
inner product can be easily deduced. For further results in this regard, including interpolation, 
we refer to [10]. 



3 Discrete Fourier analysis on the 30°— 60°— 90° triangle 

In this section we deduce a discrete Fourier analysis on the 30° -60° -90° triangle from the analysis 
on the hexagon by working with invariant functions. 



3.1 Generalized trigonometric functions 

The group A2 is generated by the reflections in the edges of the equilateral triangles inside the 
regular hexagon f2. In homogeneous coordinates, the three reflections o~i, o~2, ""3 are defined by 

tdi := -(ti,t 3 ,t 2 ), tcr 2 := — (*2,*1,*3)> ^3 '■= ~ (*3, h, h). 
Because of the relations 03 = oxa^ox = o-20~\o~2-, the group is given by 

A2 = {1, <7l, 02, a 3, <TlO~2, O-2O-1} . 

The group of isometries of the hexagonal lattice is generated by the reflections in the 
median of the equilateral triangles inside it, which can be derived from the reflection group A2 
by a rotation of 90° and is exactly the permutation group of three elements. To describe the 
elements in A\, we define the reflection —0 for any e A2 by 

t(-<j) := -t0, VtEl|. 

With this notation, the group A\ is given by 

A\ = {1, -01, -0 2 , -0 3 ,0102, 02O"l} , 

in which —01, —02, —03 serve as the three basic reflections. The group is the same as the 
permutation group £3 with three elements. 

The group G2 is exactly the composition of A2 and A\ , 

G 2 = {00* : £ A2, 0* e ^.2} = {±1> ±°"1j ±°*2) ±C3> ±0102, +02C1} • 

Let Q denote the group of A2 or A\ or G2. For a function / in homogeneous coordinates, the 
action of the group Q on / is defined by 0/(t) = f(to~), e Q. A function / is called invariant 
under Q if 0/ = / for all a e Q, and called anti-invariant under Q if 0/ = (— 1)' CT '/ for all e Q, 
where |cr| denotes the inversion of and (— 1)'°"' = 1 if = +1, +0102, +0201, and (— 1)'°"' = — 1 
if = +01, +02, +03- The following proposition is easy to verify (see [6]). 
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Proposition 3.1. Define the operators V + and V acting on /(t) by 

V ± f{t) = \ [fit) + /(t<7ia 2 ) + f(ta 2 ai) ± /(toi) ± f(ta 2 ) ± /(to 3 )] . (3.1) 

Then the operators V + and V~ are projections from the class of H-periodic functions onto 
the class of invariant, respectively anti-invariant, functions under A 2 . Furthermore, define the 
operators V* and V* acting on /(t) by 

P?f(t) = \ [/(t) + f{ta x a 2 ) + f(ta 2 ax) + /(-to x ) + f(-ta 2 ) + /(-to 3 )] . (3.2) 

D 

Then the operators V* and are projections from the class of H-periodic functions onto the 
class of invariant, respectively anti-invariant functions under A 2 . 



(-1,1,0) (0,1,-1) (-I- 1 --!) (-1,1,0) (-1,1,-5) (0.1,-1) 




(0,-1,1) (1,-1,0) (1,-1,1) (0,-1,1) (1,-1,1) (1.-1,0) 



Figure 3.1. Symmetry under A2 (left), A% (center) and G2 (right) in the physical space. The shaded 
area is the fundamental triangle of Qa under Gi- 



(-f.f.o)/ \(o,f,-f) 



(-?,<>,?)■ 



■ (f,o,-f) 



(0rS.fr 



(S--f>°) 




( 3 ' 3 > 3) 



(~ 3 >" 3 ' 3 ) 















o,-i,tN 







(!1 n _m) 
\ 3 ' 3 ' 3 ' 

(5,o,-5) 



,„ _2n ns 
v 3 ' 3 ' 3 I 



( 3 ' 3 ' 3 ) 



Figure 3.2. Symmetry under A2 (left), A% (center) and G2 (right) in the frequency space. The shaded 
area is the fundamental triangle of Qb under G2. 

For a e G 2 , the number of inversion \a\ satisfies | — a\ = \a\. The following lemma can be 
easily verified (writing down the table of era* for a e A 2 and a* e A?, if necessary). 

Lemma 3.2. Let f be a generic H-periodic function. Then 



V:V + f{t) = 1 £ (/(to) + /(-to)) = 1 2 /(to), 

n^ + /(t) = ^ 2 (/(t<r) - /(-to)) = 1 2 (-l)H (/(to) - /(-to)) , 

V+V-f(t) = 1 J] (-1)1^1 (/(to) - /(-to)) = 1 J] (/(to) - /(-to)) , 

CT6.42 cre.A* 

n~^-/(t) = ^ 2 (-i) H (/(to) + /(-to)) = 1 2 (-i)i-i/(to). 



cre.42 cret/2 
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For </>k(t) = e 3 j the action of "P + and "P on </> k are called the generalized cosine and 
generalized sine functions in [9], which are trigonometric functions given by 



C k (t) := P + k (t) 



e a 



(fcl-fc 3 )(*l-*3) 



cos k2irt2 



S k (t) := -V-fcit) 

i 6 



+ efi^k^-t,) cosk27Th + e f (fci-*B)(t3-ta) cos^vrt! 
1 



^Oi-fcaXti-ts) 



C 3 



^ sin k2irt2 



+ e f (fci-*s)(* a -*0 sin^Trtg + ef ( fc i- fc 3)fe-t 2 ) sinA;27rtl 



(3.3) 



(3.4) 



Because of the symmetry, we only need to consider these functions on the fundamental domain 
of the group A2, which is one of the equilateral triangles of the regular hexagon. These functions 
form a complete orthogonal basis on the equilateral triangle and they are the analogues of the 
cosine and sine functions on the equilateral triangle. These generalized cosine and sine functions 
are the building blocks of the discrete Fourier analysis on the equilateral triangle and subsequent 
analysis of generalized Chebyshev polynomials in [9]. 

We now define the analogue of such functions on G2. Since the fundamental domain of the 
group G2 is the 30°-60°-90° triangle, which is half of the equilateral triangle, we can relate the 
new functions to the generalized cosine and sine functions on the latter domain. There are, 
however, four families of such functions, defined as follows: 



CS k (t) 
SS k (t) 



cc k (t) := r+r + Mt) = ^ 2 (^(*) + ^(t)) = ^(c k (t) + c_ k (t)), 
sc k (t) := -/p-r + fa(t) = ^ 2 (<Mt) - 0-k*(t)) = ^(c k (t) - c_ k (t)), 

\v:v-Mt) = ~ ^ (-l)H (0 kCT (t) - 0_ k(7 (t)) = i(S k (t) - S_ k (t)), 

-r-v~Mt) = -^ 2 (-i) H (^(t) + 0- kCT (t)) = ^(s k (t) + s_ k (t)), 



where the second and the third equalities follow directly from the definition. We call these 
functions generalized trigonometric functions. As their names indicate, they are of the mixed 
type of cosine and sine functions. 

From (3.3) and (3.4), we can derive explicit formulas for these functions, which are 



CC k (t) 



SC k (t) 



CS k (t) 



SS k (t) 



1 

3L 



cos 



7r(fci-A: 3 )(ti-t 3 ) 



cos nk2t2 + cos 



COS 



7r(fci-fc 3 )(t 3 -f2) 



cos Trk2ti 



1 

3L 



sin 



7r(fci-fc 3 )(ti-t 3 ) 



cos irk2t2 + sin 



sin 



7r(fci-fc 3 )(t 3 -f 2 ) 



COS 7T/C2tl 



1 

3L 



cos 



7r(fci-A: 3 )(ti-f 3 ) 



sin7r/c2t2 + cos 



COS 



7r(fci-fc 3 )(t 3 -f2) 



sinvr/^i 



1 

3L 



sin 



7r(fci-fc 3 )fti-f 3 ) 
3 



sin irk2t2 + sin 



s n( kl 


-fc 3 )(*2 


-ti) 




3 




7r(fci- 


-k 3 )(t 2 - 


-tx) 




3 




7r(fci- 


-k 3 )(t2 


-tl) 




3 




7r(fci — 


k 3 )(t 2 - 





COS 7T/C2t3 



COS 7^2*3 



sin irk2t3 



sinvrA^s 



sin 



7r(fci-fc 3 )(t 3 -t 2 ) 



sinvr^ti 



(3.5) 



(3.6) 



(3.7) 



(3.8 
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In particular, it follows from (3.6)-(3.8) that CSk(t) = SSk(t) = whenever k contains zero 
component and SCk(t) = SSk(t) = whenever k contains equal elements. Similar formulas can 
be derived from the permutations of t\, t 2 , £3. In fact, the functions CCk and SSk are invariant 
and anti-invariant under G2, respectively, whereas the functions CSk and SCk are of the mixed 
type, with the first one invariant under A2 and anti-invariant under A 2 and the second one 
invariant under A\ and anti- invariant under A2- More precisely, these invariant properties lead 
to the following identities: 

(3.9) 
(3.10) 
(3.11) 
(3.12) 
(3.13) 

In particular, it follows from (3.6)-(3.8) that CSk(t) = SSk(t) = whenever k contains zero 
component and SCk(t) = SSk(t) = whenever k contains equal elements. Moreover, for any 
k e H', CSk(t) = SSk(t) = whenever t contains zero component and SCk(t) = SSk(t) = 
whenever t contains equal elements. 

Because of their invariant properties, we only need to consider these functions on one of the 
twelve 30°-60°-90° triangles in the hexagon fi. We shall choose the triangle as 

A := {t e R 3 H : s$ t 2 sS h «S -t 3 sc 1}. (3.14) 



CCk(ta) = 


cc k (t; 




SSk(to) = (- 


-l) |ff| SSk(t), 


0- £ Q2, 


SC k (to) = 


-sc k ( 


-to-) 


= sc k (t), 




a £ A 2 , 


CS k (tff) = 


-cs k ( 


-to) 


= (-i) H cs k 


(t), 


£ A2, 


SC k (ta) = 


-sc k ( 


-to) 


= (-i) H sc k 


(t), 


O £ A%, 


CS k (tff) = 


-cs k ( 


-to) 


= CSk(t), 




a £ Al. 



The region A and its relative position in the hexagon are depicted in Figs. 3.3 and 3.1. 




Figure 3.3. The fundamental triangles in D,a (left) and fis (right). 



When CCk, SCk, CSk, SSk are restricted to the triangle A, we only need to consider a subset 
of k £ H' as can be seen by the relations in (3.9)-(3.13). Indeed, we can restrict k to the index 
sets 

T = T cc := {k £ M t : < k 2 < fa}, T sc := {k £ M t : < k 2 < fa}, (3.15) 
T cs := {k £ H 1 : < k 2 < fa}, T ss := {kerf : < k 2 < fa}, (3.16) 

respectively, where the notation is self-explanatory; for example, r cc is the index set for CCk- 
We define an inner product on A by 

</, <?>A := £ /(t)^(t)rft = 4 £ dt 2 £ ^ f{t)W)dti. 

If fg is invariant under the group G2, then it is easy to see that (/, g)n = (/, 5)a- Consequently, 
we can deduce the orthogonality of CCk, SCk, CSk, SSk from that of 0k on ^- 
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Proposition 3.3. It holds that 



<CCk, CCj)a 

<SCk, SCj)a ; 

<CSk, CSj)a : 
(SSk, SSj)a = 



|kG 2 | 
|kG 2 | 



A 



A 



kj 



A 



kj 



|kG 2 | 

Akj 
|kG 2 | 



A 



kj 



12 "kj, 



k = 0, 

k 2 {ki - k 2 ) =0,ki> 0, 
k\ > k 2 > 0, 

k 2 = 0, 

k\ > k 2 > 0, 

h = k 2 > 0, 
k\ > k 2 > 0, 



j,keT c 

j,keT s 

J,ker c 
j,keT s 



(3.17) 

(3.18) 

(3.19) 
(3.20) 



where kG 2 = {k<r : a e G 2 } denotes the orbit o/k under G 2 . 

3.2 Discrete Fourier analysis on the 30°— 60°— 90° triangle 

Using the fact that CCk, SCk and CSk, SSk are invariant and anti-invariant under A 2 and that 
CCk, CSk and SCk, SSk are invariant and anti-invariant under A 2 , we can deduce a discrete 
orthogonality for the generalized trignometric functions. Again, we state the main result in 
terms of cubature rules. The index set for the nodes of the cubature rule is given by 

T n := {j e H : < j 2 < ji sS -j 3 < n} , 

which are located inside nA as seen by (3.14). The space of invariant functions being integrated 
exactly by the cubature rule are indexed by 

T n = := r u Mj, = {k e H T : < k 2 < h < k 3 + n}, 



-pSC 

pes 

n 

pss 



T sc u4 = {keH f 
T cs u Ml = {k e H T 
T ss u Hj t = {k£E f 



Correspondingly, we define the following subspaces of W n , 



< k 2 < k\ < k 3 + n) , 
< k 2 < ki «S k 3 + n) , 
< k 2 < k\ < k 3 + n\. 

t 



H% :=span{CC k :kE^ c }, 
K s := span{CS k :kE^}, 
It is easy to verify that 
dim?C = |C 

dim?e = irsfi 



17 s 



SCI 

n I , 



l(3Lfj-2n)(LfJ- 
r n _ 6 |, dim%" 



; span{SCk : k e T: 
span{SSk : k e r 



i) - (LfJ-«-i)(LfJ 

= \Ti c \ = dim7C = IT" 



n— 3 1 



Theorem 3.4. The following cubature is exact for all f e W^L-i 
where 

f 12, jeT° n , 



(3.21) 



(3.22) 



(n) (n)..^, | 



1, 
2, 
3, 
I 6, 



0, 

(n,0, -n), 

V 2 ' 2 ' "/ 



J 
j 
j 

otherwise, 



(interior), 
(30° -vertex), 
(60° -vertex), 
(90° -vertex), 
(boundaries). 
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(5.|>-n) 




(0,0,0) (1,0,-f) (0,0,0) (f,o,-|) 

(cs) (ss) 
Figure 3.5. The index set r„. 



Moreover, if we define the discrete inner product (f, g}& n = \ 2 ] f{ 'n)dCn) > then 

jeT„ 

A j)k A j)k 
<CCj,CC k >A,n = , n) ' = — T^y, j,ker n , 
cx kG-2 

k k 

<SCj,SC k > A ,n = { ^ 3M : = -ffi> J' ker ^ 

ci. kG 2 u;~ 

k k 

<CSj,CS k > A ,n = (n fj* - = j,ker£\ 

ci. kG 2 
k k 

<ss j5 ss k > A , n = - 7 ^- = j, k e r- 



where k = (£3 — &2, &i — ^3, ^2 — 

The formula (3.22) is derived from (2.4) by using the invariance of the functions in "H^n-i an d 
upon writing = ( k [_} tre G 2 {b(T ' * E U (Uo-eG^^ 17 : * E *^})* ^he reason that k appears 
goes back to Proposition 2.3. As the proof is similar to that in [9], we shall omit the details. 

One may note that the formulation of the result resembles a Gaussian quadrature. The 
connection will be discussed in Section 6. 
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3.3 Sturm— Liouville eigenvalue problem for the Laplace operator 

Recall the relation (2.3) between the coordinates (xi,x 2 ) and the homogeneous coordinates 
(ti,t2,ts). A quick calculation gives the expression of the Laplace operator in homogeneous 
coordinates, 



A : = 



c 2 



c 2 



1 



dx\ dx\ 2 



d_ d_ 

~ dT 2 



d_ 
dT 2 



_d_ 

dt?. 



d_ 

df 3 



d_ 



A further computation shows that </>k(t) = e s 
for kel, 



kt 



are the eigenfunctions of the Laplace operator: 



A^> k = -A k 0k, 



2V 
~9~ 



[{h - k 2 ) 2 + (k 2 - k 3 ) 2 + (h - h) 2 ] . 



(3.23) 



As a consequence, our generalized trigonometric functions are the solutions of the Sturm- 
Liouville eigenvalue problem for the Laplace operator with certain boundary conditions on the 
30°-60°-90° triangle. To be more precise, we denote the three linear segments that are the 
boundary of this triangle by B\, B 2 , B%, 



Bi 



(t e A : t 3 



1} 



Bo 



{t e A : t 2 = 0}, 



ft e A : ii 



t 2 }. 



Let denote the partial derivative in the direction of the exterior norm of A. Then 



d_ 
dn 



d_ 

dn 



B 2 



d 

dT 2 



8 

dn 



d 



Theorem 3.5. The generalized trigonometric functions CCk, SCk, CSk, SSk are the eigenfunc- 
tions of the Laplace operator, Au = —X^u, that satisfy the boundary conditions: 



CC k 

cs k 



du 
dn 
du 
dn 



B!uB 2 uB 3 



0, 



SC 



k • 



du 
dn 



BjuB 2 



0, 



u\b 3 



0, 



B 3 



0. 



u\b 1 vjB 2 = 0> SSk : "uIb^b^Bs = 0- 



Proof. Since Ak is invariant under G 2 , that is, Ak = Ak CT , Vcr e G 2 , that these functions satisfy 
Au = —XkU follows directly from their definitions. The boundary conditions can be verified 
directly via the equations (3.5), (3.6), (3.7) and (3.8). ■ 

In particular, CCk satisfies the Neumann boundary conditions and SSk satisfies the Dirichlet 
type boundary conditions. 



3.4 Product formulas for the generalized trigonometric functions 

Below we give a list of identities on the product of the generalized trigonometric functions, which 
will be needed in the following section. 

Lemma 3.6. The generalized trigonometric functions satisfy the relations, 

CCjCC k = 1 2 CC k+j(T = 1 2 CC J+k- ( 3 - 24 ) 

<reG 2 aeG 2 

CCjSC k = 1 2 SC k+j(7 = 1 2 (-l) T (SC J+lcT - SCj_ kT ), (3.25) 



o-eG 2 
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CCjCSk = — 2 CS k+j^ = ^2 2 ( CS J+kT - CSj_ kr ), (3.26) 

<reG 2 reA* 



CCjSS k = 1 2 S W = ^ 2 (-l) |r| SS j+k(T , (3.27) 

<TSG2 (TGG2 

SCjSC k = -^ 2 (-l)M(CC k+jr -CC k _ jT ) 

reA$ 

= -J2 2 (-l) M (CC j+kT -CCj_ kr ), (3.28) 

re.4* 

SCjCS k = 1 2 (-l) |T| (SS k+jT - SS k _ jT ) = 1 2 (SS j+kT - SSj_ kr ), (3.29) 

CSjCS k = -~ J] (CC k+jT - CC k _ jr ) = -1 2 ( CC J+kr - CCj_ kr ), (3.30) 

SSjSS k = 1 2 (-l)H C C k+j(T = 1 J] (-l)HcC j+kff . (3.31) 

creG2 (TGG2 

Furthermore, the following formulas hold: 

3SC 1A -i(t)CS 1) i_ 2 (t) = SS 2 ,i,_3(t), (3.32) 

[SCi, -a(t)] 2 = i[l + 2CC 1)1 _ 2 ] - [CCa^-i] 2 , (3.33) 

[CSi,i _ 2 ] 2 + [CCi.i.-J 2 = ^ [1 + 2CC 3 ,o,- 3 ] , (3.34) 

[CC10-1] 3 = — CC30-3 + 7CC10.-1 + t"CCi 1 2 + 777 + 1 _ 2 CCio,-i- (3.35) 
3d 4 o 18 2 

Proof. For (3.24)— (3.31), we only prove (3.29). Other identities can be proved similarly. By 

the definition of the generalized trigonometric functions, 

SCjCS k = 2 (- 1 ) H (^ " *-*0 x 7^ 2 " ^) 

^ 2 (-1)' T| 2 (-Ij'^'t^+J—^ + ^-Ck+j^r 



12 2 



re.4j creA* 



- <P(k.-i(rr- 1 )T ~ ^-(k-jCTT-^r] 

upon using the relation (— i)l r l+l CTT 1 \ = (— 1)H, consequently, 

SCjCS k = -^2 2 ( _1 ) kl 2 (" 1 ) |r| [^(k+ja)r + </'-(k+j ( 7)r"0(k-ja)r-^-(k-j ( T)r] 

^ 2 (-i) H (ss k+j(7 -ss k _ 



12 



proving the first equality in (3.29). Further by (3.9), 

1 2 (-l)H(SS k+jCT -SS k _ j(7 ) = l 2 (SS kCT - 1+j - SS^lj) 



~~ 12 2 (SSj+ko- - SS ko -_j) - — ^ (SSj +kcr - SSj_ kcr ), 
since SS; = SS; by (3.5). This completes the proof of (3.29). 
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We now prove the relations (3.32)-(3.35). By (3.29), 

CSi, 1 _ 2 (t)SCi, 0j -i(t) = i[(SS 2 ,i,-3(t) - SS ,-i,i(t)) + (SS_i jl>0 (t) - SS 3 ,-i,- 2 (t)) 



(SS 2i _ 2j0 (t)-SS 0j2 ,- 2 (t))l = iss 2 ,i,_ 3 (t), 



which proves (3.32). By (3.28) and (3.24), we have 



[SC 1A -i] 2 + [CCi, ,-i] 2 = -i[CC 2 ,o,-2 - 1 + 2CCi,_i >0 - 2CC 1 ,i,_ 2 ] 

+ ^[CC 2i0 ,-2 + 1 + 2CCi _i,o + 2CCi,i > _ 2 ] = g[l + 2CCi, lj _ 2 ], 
which is (3.33). Next, from (3.30) and (3.24) we deduce that 

[CSi,i_ 2 ] 2 + [CC lil5 _ 2 ] 2 = -i[CC 2 , 2 ,_ 4 - 1 + 2CC 1>li _ 2 - 2CC 3 ,o,- 3 ] 

+ ^[CC 2i2 ,_4 + 1 + 2CC lj i_ 2 + 2CC 3 ,o,-3] = g[l + 2CC 3 ,o,-3], 

which is (3.34). Finally, the identity (3.35) follows from a successive use of (3.24). The proof is 
completed. ■ 

4 Generalized Chebyshev polynomials 

In [9], the generalized cosine and sine functions Ck and Sk are shown to be polynomials under 
a change of variables, which are analogues of Chebyshev polynomials of the first and the second 
kind, respectively, in two variables. These polynomials, first studied in [6, 7], are orthogonal 
polynomials on the region bounded by the hypocycloid and they enjoy a remarkable property 
on its common zeros, which yields a rare example of the Gaussian cubature rule. 

In this section, we consider analogous polynomials related to our new generalized trigono- 
metric functions, which has a structure different from those related to Ck and Sk- 

The classical Chebyshev polynomials, T n {x), are obtained from the trigonometric functions 
cosn6 by setting x = cos 6, the lowest degree nontrivial trigonometric function. In analogy, we 
make a change of variables based on the first two nontrivial generalized cosine functions: 



x = x(t) := CCi >0 ,-i(t) = ~ (cos + cos + cos 



1 

3 

y = y(t) := CCi,i ,- 2 (t) = - (cos27rti + cos27rt 2 + cos27r(ti + t 2 )) . (4.1) 

If we change variables (t±, t 2 ) i— * (x, y), then the region A is mapped onto the region A* bounded 
by two hypocycloids, 

A* = {(x, y) : (l + 2y - 3x 2 ) (24x 3 - y 2 - 12xy - 6x - Ay - l) > 0} . (4.2) 
The curve that defined the boundary of the domain A* satisfies the following relation: 

Lemma 4.1. Let F(x,y) := (l + 2y — 3x 2 )(24x 3 — y 2 — Ylxy — 6x — 4y — 1). Then, in homogeneous 
coordinates, 

F(x, y) = 3 [SC 1A -i(t)] 2 [CSx,!,-,^)] 2 = \ [SS 2jl ,_ 3 (t)] 2 . (4.3) 



Discrete Fourier Analysis and Chebyshev Polynomials with G 2 Group 



15 




Furthermore, let J(x,y) be the Jacobian of the changing of variable (4.1); then 

Tl v 64^ 2 . . 7r(t! - i 2 ) . 7r(ti + 2t 2 ) . 7r(2*i + t 2 ) 

J(x, yj = — — sin7rri sin7ri 2 sin.7r(ii + £2) sin sin sin 

27 3 3 3 

= ^SCi,o,-i(t)CSi,i,_2(t). (4.4) 
Proof. Under the change of variables (4.1), by (3.33), (3.34) and (3.35), it follows that 
[SC li0 ,-i(t)] 2 = ^(l + 2y-3x 2 ), 

[CSi,i _ 2 (t)] 2 = 24x 3 - y 2 - 12xy -6x-4y-l, (4.5) 

from which the first equality in (4.3) follows, whereas the second one follows from (3.32). 

Taking derivatives and simplifying, we derive the formula of J(x, y) in terms of the product 
of sine functions. Furthermore, under the change of variables (4.1), it is not hard to verify that 

24rr —y — 12xy — 6x — Ay — 1 = — sin ntx sin 7ri2 sin n(ti + f 2 ), 

1 + 2y - 3z 2 = ™ sin 2 ^ ~ **> sin 2 sin 2 !^±M, 

3 3 3 3 

from which the second equality of (4.4) follows readily. ■ 

Definition 4.2. Under the change of variables (4.1), define for & 2 ^ 0, 

_ 1 _ 1 

P kuk 2 2 ( X '^) := CC fci+fc 2 ,A:2 -fei-2fc 2 ( t )) 
D |>-§/ \ SCfc 1+ fc 2+ ifc 3 fc 1 _ 2 fe 2 _i(t) 



1 1 



SCi,o,-i(t) 

CSfc 1+ fc 2 + i ) fc 2 +l,-fci-2fc 2 -2(t) 



*i,*a^.^- CSi,i_ 2 (t) 

p|.I (rr ,A ._ SS fci + fc2+2,fc 2 + l,-fc 1 -2fc 2 -3( 

n llfe2 MJ-^ ss 2> i,- 3 (t) 



We call these functions generalized Chebyshev polynomials and, in particular, call P, 2 ' 2 (2, y) 
1 1 

and P£' 2 (x,y) the first kind and the second kind, respectively. 



That these functions are indeed algebraic polynomials in x and y variables can be seen from 
the following recursive relations, which can be derived from (3.24)-(3.27). 



16 



H. Li, J. Sun and Y. Xu 



Proposition 4.3. For a, (3 = + Pj* ^ satisfy the recursion relation 

p kuk 2+ i( x >y) = Q y p kfk^ x ^y) - p kf+3M-^ x ^ ~ p kf+3M-i^y) 

p A 2+ x^y) = - p kf-, M+ ^y) ~ p kf-3,k 3+ 2^y) - ^i-i(^) M 

for k\,k2 ^ 0. Furthermore, the following symmetric relations hold, 

P%~hx,v) = K- 3 l v (x,y), PZ?„-i(x,v) = -Ct-i(^). M > 3v > 0, (4.8) 
P-Jlffav) = P~H(x,y), Pl'^ 1;V (x,y) = -pjfi^fay), v>n>Q. (4.9) 

Proof. The recursive relations (4.6) and (4.7) follow directly from (3.24) and (3.27). As for (4.8) 
and (4.9), we resort to the following identities of the trigonometric functions, 

SCn-u+i,-v,2v-fi-i(x,y) = SC( M _ 3j/ ) +v+ i ) i /)I/ _ /i _i(o;,y), 

CS[i-(u+l) + l-(v+l) + l,2v-n( x i V) = — ^(fj,-3u) + (u-l) + l,(u-l) + l,u-fi( x } V)i 
SS M _(^ + i) + 2,-( ! /+l) + l,2i/-^-l( a: ) U) = —SS(ii-3v)+(v-l) + 2,(v-l) + l,u-fi-l{ x , V)i 

^-■^-—fi+u,i',fM—2u\x,y^ = CC^+( I/ _ jU j ,v—/j,,h—2v{x., y), 

CS-^ l+ u+l,u+l,Li-2u-2(x,y) = ^^+(u-fi) + l,(iy-^) + l,^-2u-2{x,y), 

^-(fj,+ l)+u+l,u,n-2u( x ^y) = ~^{n-l) + {v-^)+l,v-n,n-2v\ X J y)i 

SS_(^+i) + !/+ 2,i/+l,^-2i/-2( a; ) y) = — ^(n-l) + {v-ii)+2,(v-Li) + l,n-2v-2{x- l J/), 

which are derived from (3.9)-(3.13). ■ 

The recursive relations (4.6) and (4.7) can be used to generate all polynomials Pj?' k recur- 
sively. The task, however, is non-trivial. Below we describe an algorithm for the recursion. Our 
starting point is 

Po,t~*{x,v) = 1, Pj'~Hx,y) = x, Po}'~H*>v) = V> 

P^p(x,y) = l, P?f*{x,y) = 6x + 2, P^p {x, y) = 6x + 3y + 1, 

p }'H x ,y) = i, Pi}'H x ,y) = 3x, P }'H x ,y) = % + 2, 

p};}(x,y) = l, P^(x,y) = 6x + 1, P^ (x, y) = 6x + 6y + 2. 

The first few cases are complicated as the right side of the (4.6) and (4.7) involve negative 
indexes, for which we need to use (4.8) and (4.9). We give these cases explicitly below 



i i 



P 2 2 ' 2 (x,y) = 6x 2 - 2x - 2y - 1, P 1 2 ' 2 [x, y) = 3xy - 6x 2 + x + 2y + 1 



i i 



P 2 2 ' 2 (x,y) = 18x 2 - 3x - 6y - 3, P 1 ^ 2 (x,y) = 18xy + 6x - 18x 2 + 6y + 3, 



i i 



P 2 2 o 2 (x, y) = 36x 2 - 6y - 3, Pf[ 2 (x, y) = 18xy + 6x + 9y + 2, 



2,0 

i i 



p>' Q * ( x , y) = 36x 2 - 6y - 3; Pf[ 2 (x, y) = 36xy + 12x + I2y + 4; 
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P30' 2 {x, y) = 36a; 3 - 18xy -9x-6y-2, 



3,0 



1 1 



P 3( p(x,y) = 108x 3 - 54xy - 27 x - 12y - 5, 



3,0 

P£i~H x , y) = 216a; 3 - 72xy - 48a; - 24y 



3,0 

P 3 3 o 5 (x, y) = 216a; 3 - 72xy - 42a; - I8y - 7; 



1 



P ,i ' 2 ( x i v) = Q y + 10 v - 72x + 36x y + 18x + 3 > 

P J'^(x,y) = 36y 2 + 36y - 216a; 3 + 108xy + 54a; + 9, 

p ofH x ,y) = 12 6a;y + !8y 2 + 36y + 54a; + 10 - 216a; 3 , 
1 1 

P Q y ( x , y) = lUxy + 36y 2 + 42y - 216a; 3 + 60a; + 11. 

The above formulas are derived from the recursive relations in the order of (2,0), (1, 1), (3,0), 
(0,2), that is, we need to deduce (3,0) before proceeding to (0,2). It should be pointed out 
that our polynomial Pq% is °f degree 3, rather than degree 2, which shows that our polynomials 

do not satisfy the property of spanjP^'^ : k\ + k^ ^ n} = II 2 . In particular, they cannot be 
ordered naturally in the graded lexicographical order. 

We shall show in the following section that our polynomials are best ordered in another 
graded order for which the order is defined by 2k\ + 3&2 = n - We have displayed the polynomials 
P k ^(x, y) for all 2k\ + 3/c2 ^ 6. In Algorithm 1 below we give an algorithm for the evaluation 

of all Pu u (x, y) with 2k\ + 3A;2 = n and n ^ 7. 



1 j_i 



The polynomials P k 2 ' 2 defined in the Definition 4.2 satisfy an orthogonality relation. Let 
us define a weight function w a ^ on the domain A*, 



(4vr 



2\a+/3 



w*A x > y) ■■= 3 2a+p ( 1 + 2 y~ 3x 2 ) (24x 3 - y 2 - 12xy - 6x - 4y - l)' 

(SC li0 ,- 1 (t))-(CS 1 , 1 ,_ 2 (t)) i 



3 

where the second equality follows from (4.5). This weight function is closely related to the 
Jacobian of the changing variables (4.1), as seen in Lemma 4.1. With respect to this weight 
function, we define 



(f,g)w aB -= c a,P f(x,y)g(x,y)w a: p(x,y)dxdy, 
J A* 



where c Qj( g := 1/ w a ^{x, y)dxdy is a normalization constant; in particular, c_i _i = 4, 

Cl 1 = c 1 1 = 18/vr 2 and ci 1 = 243/ 7r 4 . Since the change of variables (4.1) implies immedi- 

2 • 2 2 ' 2 2 ' 2 

ately that 

c a ,A f(x,y) Wa! p(x,y)dxdy=~( /(t) (SC 1 , ,i(t)) 2a+1 (CS 1)lj _ 2 (t)) 2/3+1 dt, (4.10) 
Ja* I^M Ja 

we can translate the orthogonality of CCj, SCj, CSj and SSj to that of P£ k for a, (3 = ±|. 
Indeed, from Proposition 3.3 we can deduce the following theorem. 

Theorem 4.4. For a = ±|,/3 = +\, 
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Algorithm 1. A recursive algorithm for the evaluation of P£ ^(x,y). 
Step 1 if n = 2m 

p m!o( x ^y) = 6xP m-i,o( x >y) - c P p %~2,i( x >y) - p %~2,a( x >y) - c /3 p %-3,i( x >y)> 

where cp = 2 if j3 = — |, and eg = 1 if ft = h; 
Step 2 for fc 2 from 2 — mod(n, 2) with increment 2 up to [|J — 2 do 

U _ n-3fc 2 
Kl — 2 , 



p kfkS x ^y^ = 6xP ^i,fc 2 ( a; 'y) - ^Lfe-i^'f) - p kf-2M+i( x >y^ 



Step 3 if n = 3m 

P 0,m( x > V) = % P o"m-l(>> y) - ^m-sO' V) ~ P 3°m-20> V) ~ ^m^O' V) 

+ i- p 3^-3^y) - P^-zfav), a = -h 

1 p i,m-2 ( x >y) + P l,m-1 (x, y) , a = \ ; 



if n = 3m + 1 



P 2,m-l( x > y) = GxP l,m-l( x , y) ~ P 3,m-2( x ' V) ~ P 0,m( x > V) ~ P 2,m-2( x > V) 

P l,m-l( x 'y)> a 



p ot-i^y) 



i 

2' 



0, 



a 



2 > 



if n = 3m + 2 

p i£( x >y) 



a 



3 xP t mi x >y) P 2,m-l( X ^y) P l,m—l( X ' V)i 

(6x + 1)Pq£ (x, y) - P 2 °m-i (x, y) - p i,m-l ( x >y)> a =\- 



where 




k\ = ko = 0, f -, 

h,-h i h> 0, k 2 = 0, 

, M* = 0,* 1 + fe >0, V , lS0 fe> „ 

i fei > 0, fc 2 > 0, 112 



fcl ' fc2 \i jfei > 0, fc 2 >0, klM 12' 

Proof. All four cases follow from Proposition 3.3. For a = (3 = — ^, this is immediate. For the 
other three cases, we observe that the weight function cancels the denominator in the definition 
°f P ki fc 2 P h n ( see Definition 4.2), which requires (3.32) in the case of a = @ = \. ■ 

+ i,±i 

Although the polynomials P k \ 2 are mutually orthogonal, they are not quite the usual 
orthogonal polynomials as we have seen from the recursive relations. In fact, there are only 
two such polynomials with the total degree 2, which is one less than the number of monomials 
of degree 2. As we have seen from the recursive relations, the structure of these polynomials 
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is much more complicated. To understand their structure, we study them as solutions of the 
corresponding Sturm-Liouville problem in the following section. 

5 Sturm— Liouville eigenvalue problem 
and generalized Jacobi polynomials 

Recall that our generalized trigonometric polynomials are solutions of the Sturm-Liouville 
eigenvalue problems with corresponding boundary conditions. The Laplace operator becomes 
a second-order linear differential operator in x, y variables under the change of variables (4.1). 
Using the fact that t% = —t\ — t 2 , we rewrite the change of variables (4.1) as 



1 / 27r(ti-t 2 ) , 2w(U-t 3 ) , 27r(t 3 -ti) 

X = - [COS — V 3 iJ + COS — V 3 d> + COS — — U 

1 . 

y = -(cos27rti + cos27ri2 + cos27rt3j. 
3 

A tedious but straightforward computation shows that 

(d tl -d t2 ) 2 + (d t2 -d t3 ) 2 + (d t3 -d tl ) 2 

Air 2 4tt 2 
= —[A hl (x,y)dl + 2A 1:2 (x,y)d x d y + A 2 ,2{x,y)d 2 , + 6xd x + 18yd y ] =: —£_i_i, 
y y 2 ' 2 

where we define 

A u := -6x 2 + y + 3x + 2, A X2 = A 2 i := -9xy + 18x 2 - 6y - 3. 

A 22 := -18y 2 + 108x 3 - 54xy - 27x - 9y. (5.1) 

Consequently, we can translate the Laplace equation satisfied by CCk into the equation in 

_i _i 

£ 1 1 for the polynomials P fe 2 fc ' 2 (x, y). It is easy to verify that the operator can be rewritten 
as 



£_i _i = -wi i[d x w atf s(A 11 d x + A 12 d y ) + d y uj 2> 2 (A 2 id x + A 22 d y )] 
-wi iVVi _i AV, 



where in the second line we have used 
V := (d x ,d y ) tr and A: 



tr I \ . / A H A 12 

A 2 \ A 22 



It is not difficult to verify that the matrix A is positive definite in the interior of the domain A*. 
Indeed, det A = 3F(x,y), where F is defined in Lemma 4.1, and Ai j i(x,y) = 3(x — y) + 2(1 + 
2y — 3x 2 ) is positive if x > y and it attains its minimal on the left most boundary, as seen by 
taking partial derivatives, in the rest of the domain, from which it is easy to verify that A\ \ > 

in the interior of A*. The expression of £_i _i prompts the following definition. 

2 ' 2 

Definition 5.1. For a,/3 > —1, define a second-order differential operator 

£a,/3 ■= -W- a -/3V tr W ai fsAV 

= [d x w ajj3 (A n d x + A 12 d y ) + d y w a ^(A 2l d x + A 22 d y )] . 
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5A* 



II 

A* 



The explicit formula of this differential operator is given by 

£a,/9 = -Audi - 2A 12 d x d y - A 22 d 2 y + B x d x + B 2 d y . (5.2) 

where we define 

Bi(x, y) = Tlx + Ylax + 18/3x + 6a + 3, 

B 2 {x, y) = \%x + 36ax + 18/3 + 45y + 36/3y + 18ay + 9. 

Theorem 5.2. Let a,/3 > —1. Then, the differential operator 

£a,P = -lO-a.-^V'V^AV 

is self-adjoint and positive definite with respect to the inner product (•, -) Wa „. 
Proof. By Green's formula, 

JJ f£a,pgw<x,pdxdy = JJ fV tr w a)/3 AVgdxdy = - ^ (V f) tr A(V g)w a)/3 dxdy 
A* A* A* 

<j> w a ,pf [{Aiid x g + A 12 d y g)dy - f{A 22 d y g + ^21^.5)^] 
(Vf) tr A(Vg)w a ,pdxdy 

+ j> w a ,pf [{dx9){M\dy - A 2X dx) - {d y g){A 22 dx - A 12 dy)] , 

5A* 

where dA* denotes the boundary of the triangle. Recall that 5A* is defined by F(x,y) = 0, 
where F is defined in Lemma 4.1. It follows then 

dF dF 

dF = F x dx + F 2 dy = 0, where F x = — , F 2 = — . (5.3) 

ox oy 

On the other other hand, a quick computation shows that 

FiA u + F 2 A 21 = -6(5x + l)F(x, y) = 0, (5.4) 
F X A X2 + F 2 A 22 = — 6(3j/ + 2x + l)F(x, y) = (5.5) 

on dA* . Solving (5.3) and (5.4) shows that Andy — A 2 \dx = 0, whereas solving (5.3) and (5.5) 
shows that A 22 dx — A\ 2 dy = on dA* . Consequently, the integral over dA* is zero and we 
conclude that 

- f£a,pgwa,pdxdy = JJ(V ff 'A(V g)w atf3 dxdy = - ^ gC^pfwa^dxdy, 
A* A* A* 

which shows that is self-adjoint and positive definite. ■ 

We consider polynomial solutions for the eigenvalue problem C a ^u = Xu. Differential opera- 
tors in the form of (5.2) have long been associated with orthogonal polynomials of two variables 
(see, for example, [8, 16]). However, in most of the studies, the coefficients A{j are chosen to 
be polynomials of degree 2, which is necessary if, for each positive integer n, the solution of the 
eigenvalue problem is required to consist of n + 1 linearly independent polynomials of degree 
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n, since such choices ensure that the differential operator preserves the degree of polynomials. 
In our case, however, the coefficient ^2,2 in (5.1) is of degree 3, which causes a number of 
complications. In particular, our differential operator does not preserve the polynomial degree; 
in other words, it does not map Ii 2 n to nl, the space of polynomials of degree at most n in two 
variables. 

Definition 5.3. For k\,k 2 > 0, the m-degree of the monomial x kl y k2 is defined as : = 
2k\ + 3k 2 . A polynomial p in two variables is said to have m-degree n if one monomial in p has 
m-degree of exactly n and all other monomials in p have m-degree at most n. For n e No, let II* 
denote the space of polynomials of m-degree at most n; that is, 

IT; := span{x fcl /' 2 :0^k 1 ,k 2 ; 2ki + 3k 2 ^ n}. 

The dimension of the space II* is the same as that of H 1 ^, by (3.21), 

dim U* = \ (3Lf J - 2n) Qf J + l) - (Lf J - n - l) Qf J + l) . (5.6) 

Here is a list of the dimension for small n: 



n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


dim II* 


1 


2 


3 


4 


5 


7 


8 


10 


12 


14 


16 


19 



The name m-degree is coined in [13] after the marks, or co-marks, in the root system for 
the simple compact Lie group, where the case of the group G 2 is used as an example. For 
polynomials graded by the m-degree, we introduce an ordering among monomials. 

Definition 5.4. For any k,j 6 Njj, we define an order < by k < j if 2j\ + 3j 2 > 2k\ + 3k 2 
or 2(k\ — ji) = 3(j2 — k 2 ) > 0, and k<jifk<joik = j. We call < the *-order. If 
Ck\,k2 x ^ x y^ 2 with c mtn 0, we call c m ^ n x m y n the leading term of p in the 

(ki,k2)<(m,n) 

*-order. 

For m, n > 0, define 

n m,n = span{x J y k : (j, k) < (m,n)}. 
It is easy to see that II* = II* . 2n . . 2 „ , ■ The *-order is well-defined. The following lemma 

L 3 J ' L 3 J ^ 

justifies our definitions. 

Lemma 5.5. For m, n > ; the operator C a ^ maps n^ n onto I^m,n- 

Proof. We apply the operator C a ,i3 on the monomial x^y k . The result is 

£ a ,px j V k = -A n d 2 x x j y k - A 22 d 2 y x j y k - 2A 12 d x d y x j y k + B 1 d x x j y k + B 2 d y x j y k 
= [6(j 2 + 3k 2 + 3jk) + 3(5 + 4a + 6(3)j + 3(9 + 6a + 120) Jfe] x j y k 

- 108k(k - l)x j+3 y k - 2 - j(j - l)x j - 2 y k+1 + I8k(3k - 2 - 2j + 2a)x j+1 y k - 1 
+ 3j(~j + 2 + 4k + 2a)x j - 1 y k + 9k(k + 2f3)x j y k - 1 

- 2j(j - l)x j - 2 y k + 27k(k - l)x j+1 y k - 2 + 6jkx j - 1 y k - 1 . 

Introducing the notation 

T = {(0, 0), (0, 1), (1, 0), (1, 1), (2, 1), (3, 2), (4, 2), (4, 3), (5, 3)} , 
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we write the expression as 

C a ,^y k = 2 a^- 2 ^y k+ ^, (5.7) 



where 



oj* = 6(j 2 + 3/c 2 + 3j/c) + 3(5 + 4a + 6/3) j + 3(9 + 6a + 12/3) A;, 
oj* = -108Jfe(fc - 1), oj'5 = -j'(j - 1), oi'* = 18fc(3fc - 2 - 2j + 2a), 



•!>* = 3j{-j + 2 + 4k + 2a), afj = 9fc(fc + 2/3), 

g = -2j(i-l), ag = 27fc(fc - 1), 4 k 3 = 6jk. 



From this computation, it follows readily that C a p maps n^ n into n^ n . Furthermore, with 
respect to the *-order, it is easy to see that a™£ l x m y n is the leading term of C a> p by (5.7), which 
shows that £ Qj/ g maps 11^ n onto 11^ n . ■ 



The identity (5.7) also shows that C a ,p has a complete set of eigenfunctions in II,; 



* 



Theorem 5.6. For a,j3 > —1/2 and k\,k 2 5= 0, i/iere exists a polynomial P^fk 2 E ^fci ife 
the leading term x kl y k2 such that 

where 

K[% ■= + 5 + 4a + 6/3) + h 2 (k 2 + 1 + 2/3). (5.9) 

Furthermore, if we require all the polynomials are orthogonal to each other with respect to the 
inner product (•, •) U)q g , then P^fk 2 ^ s un ^Q ue ^y determined by its leading term in the *-order. 

Proof. We first apply the Gram-Schmidt orthogonality process to monomials {x fcl y fc2 } in the 
*-order, which uniquely determines a complete system of orthogonal polynomials with leading 
term x kl y ki with respect to (-, -} Wa g ; that is, Pqq f(x, y) = 1 and 

rtfto&v) = * kl v k2 ~ E < fX , ^ Wa ' P ^ ix ^ (0 ' 0) < {klM) - 

The Gram-Schmidt orthogonality and Lemma 5.5 show that 

£a,pPZS 2 (*,y) e span{P«'J 2 (x,y) : (j u j 2 ) < {h,k 2 )} = H* klM . (5.10) 
Evidently, C^Pqq = = Oq 'qPqq j . We apply induction. Assume that 

= °°i J 2 > 0i , J2) < (fci , k 2 ) . 
It then follows from Theorem 5.2 and the orthogonality of -Pj^' j. that 

so that, as a consequence of (5.10), 
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Comparing the leading term of the above identity, we obtain from (5.7) that 



Jfci>&2_fci,fc2 _ ™fci,fca 



which gives c = Oq 1 ^ 2 . Ultimately, this inductive process shows that 

L a,^ kl> k 2 - A fci,Jfe2 i fci,fc 2 ' Wlt±1 A fci,fc 2 - °0,0 • 
As shown in the proof of Lemma 5.5, 

A £f*a = 6 ( fc i + 3/c 2 + 3hk 2 ) + 3(5 + 4a + 6/9) fci + 3(9 + 6a + 12/3)& 2 

= -(2k! + 3fc 2 )((2fci + 3fc 2 ) + 4 + 4a + 6/3) + -fc 2 (fc 2 + 2 + 2/3) + 3fci, 

which is (5.9) since \k\* = 2k± + 3/c 2 by definition. 

Moreover, suppose Pj? k (x,y) e LT^ fca is another polynomial with the leading term x kl y k2 
such that 

C a,pPkS 2 (a?, y) = A-P£'J 2 (x, y), 

= °' V ^ E s P an {^V 2 : (iii J2) < (fei, ^2)} • 

Using the same argument that determines c in (5.11), we see that A = A]j^ ^ = floo^- Moreover, 
it is easy to see that 

p kS 2 - Km e s P an i xjl y j2 ■ &>h) < (**> ^) > 

This finally leads to P^fk 2 ~ Pkfk 2 = ^> which shows that P£ ^ is uniquely determined by its 
leading term in the *-order and the orthogonality (P£ k , x^ 1 y^ 2 ) Wa ^ = for all (j\ , j 2 ) < ( fci, & 2 ) . 
This completes the proof. ■ 

Let Pfa fa be orthogonal to each other with respect to the inner product (•, -} Wa /3 . The first 
few polynomials and the eigenvalues can be readily checked to be 

'/ „,\ _ 1 \ a 'P 



P^(x,y) = l, A o 7 = 0; 
Ptf(x,y) = x + 7+ \ + J^ 6/r Agf = 3(7 + 4a + 6/3); 

a3 3(1 + 2a) 5 + 5a + 11/3 + 2a/3 + 6/3 2 + 4a 2 

°* {X,y > ~ V+ 4 + a + 3/3 X+ (4 + a + 3/3)(5 + 2a + 4/3) 
Aof = 9(5 + 2a + 4/3); 
P^fziil-rr 2 2 ^ „ ■ 4(a + l)(2a-l) _ 



3(3 + 2a) * (3 + 2a)(4a + ll + 6/3) 
-105 - 86a - 120/3 - 36,9 2 - 48,9a + 8a 2 + 24a 3 
+ 3(3 + 2a)(4a + 11 + 6/3)(4a + 6/3 + 9) ' 

A 2 ,o = 6(9 + 4/3 + 6a), 
afB 3(2a-l) 2 6/3a + lla + 15/3 + 2a 2 + 27 

x ' 1 {X,y > ~ Xy+ 5 + a + 3/3 X + (5 + a + 3/3)(4a + 6/3 + 13) y 

1 19 + 229/3 + 40a 3 + 36/3 3 + 1 1 la + 80/3a 2 + 156/3 2 + 132/3a + 140a 2 + 36/3 2 a 
(2a + 7 + 4/3)(4a + 6/3 + 13)(5 + a + 3/3) 
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r»„ :=\(m-2p + 3q,n+p-2q)eZ 2 : < g < _l, 1 < p < [m^j I 



8a 3 + 6a 2 + 4/3a 2 + 12/3 2 a + 9a + 28/3a + 70 + 93/3 + 30/3 2 
+ (2a + 7 + 4/3)(4a + 6/3 + 13)(5 + a + 3/3) ' 

A^f = 6(14 + 9/3 + 5a). 

For each Pm,n, (5.10) shows that the jC a ,pPm,n involves only Pf^j 2 with (ji,j 2 ) in 
r m ,„ := {(ii, j 2 ) 6 N 2 : (ji,j 2 ) < (m,n)} . 

This set of dependence of the polynomial solution is determined by the *-ordering. Indeed, it is 
easy to see that 

r m,n : = {{m-2p + 3q,n + p-2q)) e Z 2 : s; q < [2^J,0 < 2m + 3n} , 

For p, q as in but not both 0, we have that for a, (3 Ss — ^, 

A fci 2 " A ^2 P+ 3 9 , fc2+P -2 9 = 3(2fci - 2p + 3q + 2a + l)p + 9(2fc 2 + p - 2q + 2/3 + l)g > 0, 

which shows that ^ k f k2 ^jfj 2 f° r an y {31,32) E rjj" fca . This implies that polynomial solutions 
of the same m-degree below to different eigenvalues. Moreover, if \ k f k2 = ^jfj 2 f° r {hi 32) < 
(ki,k 2 ), then {31,32) e T kiM . 

In the case of (a, /3) = (— |, — i), our polynomials i^'^ 2 agree with the generalized Chebyshev 
polynomial that we defined in the last section. For the other three cases of (a, /3) = (ijj+s), 
this requires proof. Let us denote the Chebyshev polynomials temporarily by Q k f k , (a,/3) = 
(— \,— |). It is not hard to see, from Algorithm 1, that the leading term of Q k f k2 is cx kl y k2 

with certain c > 0, which implies that spanjQ^' j 2 (x, y) : {31,32) ^ (^1,^2)} = k . Thus, we 
can write 

(j 1 j2)<(fci,fc 2 ) 

On the other hand, by the orthogonality and the self-adjointness of C a ^, for any {h,l 2 ) < 

{h,k 2 ), 



Q a,P a,p y h,l 2 Q a,fi \ r> 

(jiJa)<(ii,ia) / Wa/3 

As a result, we deduce from (5.12) that 

£*,f3Qtfk 2 { x >y) = KffoQkfM&y}- 

Consequently, up to a constant multiple, we see that Q k f k2 coincides with the Jacobi polyno- 
mials. 

Corollary 5.7. The Chebyshev polynomials defined in Definition 4.2 satisfy the equation (5.8). 

In particular, this shows that the Chebyshev polynomials are elements in nfL and they are 
determined, as eigenfunctions of £ a ,p, uniquely by the leading term in the *-order. 
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6 Cubature rules for polynomials 

In the case of the equilateral triangle, the cubature rules for the trigonometric functions are 
transformed into cubature rules of high quality for polynomials on the region bounded by the 
Steiner's hypocycloid. In this section we discuss analogous results for the cubature rules in the 
Section 3. To put the results in perspective, let us first recall the relevant background. 

Let w be a nonnegative weight function defined on a compact set Q in M 2 . A cubature rule 
of degree 2n — 1 for the integral with respect to w is a sum of point evaluations that satisfies 

f r 

f{x)w{x)dx = V Xjf(xj), Xj e R 

for every / e n^—i- It is well-known that a cubature rule of degree 2n — 1 exists only if N ^ 
dimll^_ 1 = n(n + l)/2. A cubature that attains such a lower bound is called Gaussian. Unlike 
one variable, the Gaussian cubature rule exists rarely and it exists if and only if the corresponding 
orthogonal polynomials of degree n, all n + 1 linearly independent ones, have n{n + l)/2 real 
distinct common zeros. We refer to [4, 15] for these results and further discussions. At the 
moment there are only two regions with weight functions that admit the Gaussian cubature 
rule. One is the region bounded by the Steiner's hypocycloid and the Gaussian cubature rule is 
obtained by transformation from one cubature rule for trigonometric functions on the equilateral 
triangle. 

6.1 Gaussian cubature rule of ra-degree 

We first consider the case of wi 1 , which turns out to admit the Gaussian cubature rule in the 

2 ' 2 

sense of m-degree. 

Theorem 6.1. For wi 1 on A* . the cubature rule 

2 ' 2 

ci 1 JJ f(x, y)whk ( x , y )dxdy = * 2 J] |SS 2 a,-3(^)| 2 /(x(^), y(^s)), (6-1) 

A* J e ^n + 5 

is exact for all polynomials f e P£ n -i ■ 

Proof. Using (4.10) with a = (3 = \ and (3.32), we see that 

ci 1 f f{x, y)wi 1 (x, y)dxdy = f f(x(t), y(t)) [SS 2 ,i,_ 3 (t)] 2 dt. (6.2) 
2 ' 2 Ja* 2 ' 2 l A l Ja 

By (4.3) and (4.5), [SS 2 ,i,- 3 (t)] 2 has m-degree 10, so that /(x(t),y(t)) [SS 2i i,_ 3 (t)] 2 6 Uf n+9 if 
/ e n^ n _i- Since SS2,i,-3(t) vanishes on the boundary of A, applying the cubature rule (3.22) 
of degree 2n + 9 to the right hand side of (6.2) gives the stated result. ■ 

What makes this result interesting is the fact that, by (3.21), 

I 1 = 1^+5! = |r„ c _i| = dimll*_ 1 , 

which shows that the cubature rule (6.1) resembles the Gaussian cubature rule under the m- 
degree. Furthermore, it turns out that it is again characterized by the common zeros of ortho- 
gonal polynomials. Let Y° be the image of : j 6 T^ s } under the mapping t >—>■ x, 

1 1 

which is the set of nodes for (6.1). Then all polynomials P£'\ with m-degree n vanish on Y°. 
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Theorem 6.2. The set Y° +5 is the variety of the polynomial ideal \Pfa% 2 (%) ■ 2&i + 3&2 = n). 

i i 

Proof. By the definition of P£'% , it suffices to show that 

SS k (^§) =0 for jeT, kef and k 1 -k 3 = n + 5. (6.3) 



Directly form its definition, 



SSk(^) = 3 



sin ^(fci"fc3)0i- J3 ) in ^km + sin Aki-k 3 )(h-h) in ^ 

3(n+5) n+5 3(n+5) n+5 

„;„ 7r(fci-fc 3 )03-j 2 ) • 7rfc 2 ji 
blH 3(n+5) bH1 n+5 



- sin 
3L 



M sir. 7rfc2 J' 2 _1_ sin ^O'a-Jl) „;„ ^2,73 , „;„ ^(j3-j2) ; 7rfc 2 jl 
3 bH1 n+5 ~r blH 3 n+5 ~r blIi 3 5,111 n+5 



Since j\ = j'2 = j'3 (mod3), we conclude then SSk(^L^) = 0. The proof is completed. 



In [13] , the existence of the Gaussian cubature rule in the sense of m-degree and the connection 
to orthogonal polynomials were established in the context of compact simple Lie groups. The 
case of the group G2 was used as an example, where a numerical example was given. The 
domain A* and the one in [13] differ by an affine change of variables. 

Our results give explicit nodes and weights of the cubature rule and provide further expla- 
nation for the result. 

(-3.1) y /(i.i) \ / (i.i) 



(44 




Chebyshev-Guass 

(-I.i) \ / ' 1 




(44) 



Chebyshev-Guass-Lobatto 

(- 5 ,i)\ /(i.i) 




(4.4) 




(44) 



Chebyshev-Guass-Radau I Chebyshev-Guass-Radau II 

Figure 6.1. The cubature nodes on the region A*. 



6.2 Gauss— Lobatto cubature and Chebyshev polynomials of the first kind 

In the case of w_i _i, the change of variables t >—* x shows that (3.22) leads to a cubature of 

2 • 2 

m-degree 2n — 1 based on the nodes of Y n . 
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Theorem 6.3. For the weight function w_i _i on A* the cubature rule 



c_i _l 
2 • 2 



A* 

/joZefe /or / e II| n _ 1 . 



^f(x,y)w_i _i{x,y)dxdy = ^ J] u;i n) f(x(±), y^)) 

A * 

jeT n 



(6.4) 



The set Y„ includes points on the boundary of A*, hence, the cubature rule in (6.4) is an 

analogue of the Gauss-Lobatto type cubature for w_i _i on A*. The number of nodes of this 

2 > 2 

cubature is dim II*, instead of dim 11*^ . In this case, the corresponding orthogonal polynomials 

_i _i 

are the generalized Chebyshev polynomials of the first kind, Tk ly k 2 {x,y) := P n 2 ' 2 {x,y). The 
polynomials in {T a : |a|* = n} do not have enough common zeros in general. In fact, the two 
orthogonal polynomials of m-degree 6, 

T 3 , (x, y) = 36x 3 - 18xy - 9x - 6y - 2, 
T 2 , 2 (x, y) = 6y 2 + lOy - 72x 3 + 36xy + 18x + 3. 



only have three common zeros on A*, 

(x,y) 



#7Cos(^ 



i arccos _ihg_; 



V7+1 



M = 0,l,2, 



whereas dim IT*, = 5. For cubature rules in the ordinary sense, that is, with 11^ in place of II*, 
the nodes of a cubature rule of degree 2n — 1 with dim H"^ nodes must be the variety of a poly- 
nomial ideal generated by dimIX* +1 linearly independent polynomials of degree n + 1, and these 
polynomials are necessarily quasi-orthogonal in the sense that they are orthogonal to all polyno- 
mials of degree n — 2 [19]. Our next theorem shows that this characterization of such a cubature 
carries over to the case of m-degree. 



Theorem 6.4. Denote a* = (cti — 1,02); &i > a 2 , a nd a* 

is the variety of the polynomial ideal 

(T a (x) - T a *(x) : |a|* = n + l>. 



(«i, cti — 1) if ax = ct2- ThenY n 



(6.5) 



Furthermore, the polynomial T a (x)—T a *(x) is of m-degree n+1 and orthogonal to all polynomials 
in II* _ 2 with respect to w_ 1 _ 1 . 



Proof. A direct computation shows that, for any ke I with k\ 



n 



CCfc 1 _i i fc 2; fc 3+ i(t) — CC k (t) 



1 

3L 



cos 



7r(n-l)(ti-t 3 ) 



; I ; 1 ir(n— l)(t'>— U) 1 , 

— COS !TK2t2 + COS — COS 7TK2I3 



COS 



ir(n-l)(t 3 -t 2 ) 



COS Tlk2t 



COS 



7r(n+l)(ti-t 3 ) 



cos 7rA;2t2 



COS 



7T(n+l)(t 2 -h) 



., cos irk2t3 + cos 7r ("+ 1 K t « ^ cos 7T ) (2 t 1 

sin — ^ — — sin v 3 ' cos 7rfc2i2 + sin — L -| — — sm 3 ; cos -KK2H 
+ sm — — — sm v 3 ' cos 7r/c2ti , 
where we have used the definition of CC^ for the first equality sign. Hence, for any j e T n , 



CC 



fcl — 1,^2,^3 + 1 



cc k 



sm 



sm 



3n 



COS 
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sin u , ^ ; sin 



3n 



COS 



bill g bill 3n COb n 



where the last equality sign uses the fact j\ = j'2 = j'3 (mod 3). With ot\ = k\ + ^2, this shows 
that T a — T a * vanishes on Y n . Finally, we note that = | ck | ^ — 2 or |a|* — 1, so that T a * is 
a Chebyshev polynomial of degree at least n — 1 and T a — T a * is orthogonal to all polynomials 

in n*_ 2 . ■ 



6.3 Gauss— Radau cubature and Chebyshev polynomials of mixed kinds 

Under the change of variables t >—* x defined in (4.1), we can also transform (3.22) into cubature 

rules with respect to w_i 1 and w_i 1, which have nodes on part of the boundary and are 

2 ' 2 2 ' 2 

analogue of Gauss-Radau cubature rule. They are associated with Chebyshev polynomials of 
the mixed types. We state the result without proof. 

Theorem 6.5. The following cubature rules hold, 



/Hn^)*^))* v/en* n _!,(6.6) 



M£s)>v{i&3))> V/e 11^.(6.7) 

Since by (4.5), SCi^-i and CSi^-2 vanish on part of the boundary of A, the summation is 
not over the entire T n +2 or T n +3 but over a subset that exclude points on the respective boun- 
dary. Let Y^+! and Y£+ 3 denote the set of nodes for the above two cubature rules, respectively. 

Theorem 6.6. Y^_ 2 ^ s ^ e variety of the polynomial ideal 

(Pa*'*(x) : H* =n>. (6.8) 
And is the variety of the polynomial ideal 

<pj'~^(x) -pjH(z) : |a|* = n+ 1>. (6.9) 

It is of some interests to notice that, in terms of the number of nodes vs the degree, (6.6) is 
an analogue of the Gauss cubature rule in m-degree. 
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